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Abstract. In this paper we study the (asymptotic and exponential) stability 
of the m-fold circle as a solution of the p-curve shortening flow (p > 1 an 
integer) . 
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1. Introduction 

Let 

a; : X [0,r) — > 

be a family of smooth immersions of S^, the unit circle, into R^. In this paper we 
will say that x satisfies the p-curve shortening flow, p> \, ii x satisfies 

dx 1 

where k is the curvature of the embedding and N is the normal vector pointing 
outwards the region bounded by x{-,t). 

Much is known about this family of flows. To select a few among many beautiful 
and fundamental works on the subject, we must mention the works of Gage and 
Hamilton ([7]), and of Ben Andrews ([1[3]). 

In this paper we will be concerned with the stability of m-fold circles as solutions 
to the p-curves shortening fiow. As it is well known there are small perturbations 
of the 2-fold circle that do not behave asymptotically as a shrinking 2-fold circle. 
However, very recently, Wang in [11] showed the asymptotic stability of m-fold 
circles under certain small ^^-periodic perturbations as solutions to the curve 
shortening flow (i.e., for the case when p = 1). In this note we will extend the 
work of Wang in two ways: we will show asymptotic stability results for the m- 
fold circle as a solution to the p-curve shortening flow for p any positive integer, 
and we will provide sharp stabilization estimates for the curvature of solutions to 
the p-curve shortening flow that are appropriate small perturbations of an m-fold 
circle. Other interesting works, besides Wang's, regarding stability of solutions to 
the curve shortening flow are the by now classical papers of Abresch and Langer 
([I]), and of Epstein and Weinstein ([5]). 
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To study the stability of m-fold circles as solutions to (HI we will consider the 
Boundary Value Problem 



(2) 



ip- 



on 



1) kP-'^k^g 



ifcf) in [0,f]x(0,T) 



with periodic boundary conditions, A > y ^-nd tp ^ strictly positive function. 
As it is ^ has no immediate geometric interpretation. However, when A is an 
appropriate rational number, ([2]) is the evolution equation of the curvature of a 
curve being deformed via ^ under the assumption that the curve being deformed 
satisfies certain symmetries; more precisely, when A = — , the study of equation ([2]) 
is equivalent to the study of ([T]) when the initial data is a perturbation of an m-fold 
circle under a 2^-periodic perturbation. 

The method we will use to prove our stability results was introduced in [4 (in- 
spired by [9* ) to study the blow-up behavior of certain nonlinear parabolic equations 
with periodic boundary conditions, but as the reader will notice, it can be also used 
to study the stability of certain blow-up profiles, and the regularity of solutions to 
([2]) (as a byproduct of the method we will employ, it can be shown that, under 
certain conditions on the initial data, solutions to ^ are analytic). So we hope 
that the reader may find the method used in this paper of independent interest. 

The organization of this paper is as follows: in Section [2] we present our main 
result, its application to the stability problem of m-fold circles, and its proof; in 
Section [3] we discuss the exponential stability of the constant steady solution of the 
normalized version of 

2. Main Result 

Our results on the behavior of the p-curve shortening fiow will follow as a conse- 
quence of a result on the behavior of solutions to ([2|) , that we will promptly describe; 
but before we state our main result, let us set some definitions and notation. Given 
/ G ([O, x])' write its Fourier expansion as, 

^/(n)e*^"^ 
and define the family of seminorms, 



l/ll 



max < sup \n\ 



Re (/(n)) 


,sup \nf 











As is customary, we define ([O, ^]) as the space of functions with continuous 
derivatives of order Z, equipped with the norm 



sup 



'SG 0, 



The main result of this paper is the following theorem: 
Theorem 2.1. Let A > y There exists a constant Cp^x > such that if 



(3) 



4^(0) de>cp,x\m 
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then a solution to (0) with initial condition V' is analytic (in 9) and satisfies 

(4) k{9,t)~ki0,t) ^ ^ , <i?i,p.A(T-t)(^'-'^)^, 

where < T < oo is the blow-up time of the solution to with initial condition 

v.. 

Let us remark that the fact that a solution to ^ with initial condition V' > 
blows up in finite time is a consequence of the Maximum Principle for Parabolic 
Equations. Also, the reader should notice that if Cp^\ > is large enough then there 
is no need to assume that is strictly positive, since it would be a consequence of 
((31): this is why even though we are assuming the positivity of ip, this assumption 
does not appear explicitly in the statement of the theorem. 

Next we give the promised geometric corollaries of our main result. 

2.1. A stability result for m-fold circles. Theorem 12.11 has as a corollary an 
asymptotic (nonlinear) stability result for small perturbations of m-fold circles. In- 
deed, if an rn-fold circle is perturbed by a ^^p^-periodic function, then its curvature 
function satisfies equation ^ with X — n and m mutually primes. Hence, if 



n j p - 



m \j p 

and if the perturbation of the m-fold circle by the ^^-periodic function is such that 
its curvature function satisfies the hypothesis of Theorem 12.11 then it will shrink 
asymptotically as an m-fold circle under the p-curve shortening flow. 

Let us explain with more care. Adopting a similar notation as in [111, let 7^ be 
the m-fold circle, i.e., 

7„ (6») = (cos 6*, sin 6*) , 6* e [0, 27rm] , 

and let (/? be a smooth function of the normal angle {9) of 7,„ of period with 
— > Consider the curve 

m \ P 

where N is the outward unit normal to 7^. It is assumed that \5\ > is small 
enough so that the curvature of 795 is strictly positive. Then we have the following 
result: 

Theorem 2.2. // \6\ > is small enough, then the solution to the p-curve short- 
ening flow shrinks to a point asymptotically like an m-fold circle. 

How small is \S\ in the previous theorem is dictated by the fact that the curvature 
of 7o5 must satisfy the hypothesis of Theorem 12.11 The proof of this result is 
inmediate from Theorem 12. 11 and it extends the results of Wang ([11 J, at least in 

the case when ^ > \J^^- Notice also that the condition ^ > \J^^ translates 

to ^ > \/3 in the case of the curve shortening flow which is obviously weaker than 
the Abresch-Langer condition ^ > ^/2 (which is covered in Wang's work). 

We can be more quantitative in describing how a solution to the normalized 
p-curve shortening flow with initial condition 705 approaches the solution given by 
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an m-fold circle. To this end we consider the foUowing normahzed version of ([2|) 



(5) 




pk^+p -pk in [0,^] X (0,oo) 



P+i I ^' 



with periodic boundary conditions. This normahzed equation is obtained from ([2]) 
by the rescaUng and change of time parameter given by 



k{e,t) = 

So we have: 



p+l 



(T-t)— ki9,t), 



1 



p+l 



log ( 1 - ^ 



Corollary 2.1. Let X 



n and m relatively primes, and A > Let 705 be 



2tt 



periodic perturbation of the m-fold circle. There exists a Cp^\ > such that if 



\6\ is small enough so that the curvature 0/705 satisfies 



(6) 



A 

2^ 



1/; (0) de 



2^ 



2 ' 



then the curvature k(-,T) of the solution to the normalized p-curve shortening flow 

1 

with initial condition (^^^^ ''^ I06 satisfies 



k{0,T) 



2^ 



k {9, t) dO 



< Ai^p^x exp {-13 {p, I) t) , 



Co, 



where 



/? (p, A) = A^ - 



P 



P+l = X^p — p — 1, 



and Ai,p,\ > is a constant, which only depends on the initial condition, and also 
on p and A. Furthermore, the deformation towards a circle is through analytic 
curves. 

The stabilization rate given by the previous corollary is optimal in the sense that 

p + 2' 



/3(A,p)= A^- 



P 



p+l = Xp — p — 1, 



is the smallest eigenvalue of the linearization around w = 1 of the nonlinear operator 



(7) 



pu 



pu 



pu, 



0, 



27r 



which corresponds to the elliptic part of the parabolic equation obtained by normal- 
izing ^ as is described in the introduction, so we cannot expect better stabilization 
rates for the derivatives of k. In Section [31 the last one of this paper, we also show 
that in the case of the normalized flow fc — )■ 1 at the rate predicted by the lineariza- 
tion of (O: standard methods predict that this rate is at least e~*^'^^^'^^~'^^'^, for any 
e > 0; the nonobvious part is then to dispose of the e. The reader is advised to 
consult the interesting work of Wiegner on the subject (see |12])- 
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2.2. Proof of Theorem 12.11 The proof of Theorem 12.11 will be given through a 
series of lemmas, but the strategy we will follow can be described very succinctly: 
we shall show that the Fourier coefficients of the solution to ([2]) minus its average 
decay uniformly to in time at the appropiate rate, when the blow-up time is 
approached. Our main tool is the analysis of the infinite dimensional ODE system 
satisfied by these Fourier coefficients. 

To start with the proof of Theorem 12. !( let us introduce some notation to make 
our writing a bit easier. 

u*™ (gi, (72, • ■ ■ ,qm,t) ^ it iqi.t) u {q2,t) .. .u {qm,t) , 

H{p,qi,q2) = - - (p- 1) X^qiq2 - X^qf, 
P 

and 

q= {qi,q2,---,qp+i,qp+2)- 

Let Z be a finite subset of the integers which contains (i.e, £ Z), and which is 
symmetric around (i.e., if n e Z then — n G Z). Let 

Bn^ {{bi,..., bp+2) e Z^^^ : bp+2 = 71-61-62 6p+i } , 

and define the set 

An = {b e Bn ■ there are l<i<j<p + 2 such that bi ^ and bj ^ 0} . 

Consider the following finite dimensional approximation of ([2]) (obtained from ([2|) 
after formally taking Fourier transform and then restricting the infinite dimensional 
ODE system only to those Fourier wave numbers contained in Z), 
(8) 

f ikz{0,t) = fkz{0,tr+' + 

T,qeAonZp+^ H (P, <?!, '72) k*^''^"^^ (q, t) , 

ikz{n,t) = [p±l-\^n^ykz{0,tf+^kzM + 

Eqeyt„n2p+2 H {p, gi, (J2) k*J''^^^ (q, t) , if n ^ 0, n € Z, 
with initial condition 

(9) i^z (n) = V' (") , if n e Z. 

Notice, and this will be important but not explicitly mentioned in our arguments, 
that the symmetry of Z guarantees that 

^ H {p,qi,q2)k*J^^^^ {qi, . . . , qp+2,t) 
qeAonZp+^ 

is real valued. 

Our first lemma gives an interesting estimate on the behavior of solutions to ([8|) . 

Lemma 2.1 (Trapping Lemma). There exists a constant Cp^x > independent of 
the choice of Z such that if the initial datum ip satisfies @) then there exist a 7 > 
that depends on -0 such that the solution to (0) satisifies 



kz (n, t) 



Cp>(0)e-Tl"l* 
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Proof. Following [4], first we fix a set Z as described above. Define the function 

V (n, t) as 

V {n, t) := tp {n,t) kz {n, t) , 

where, 

with 7 small, and we let 

$ gi,...,(7p+2,< = —f Tn :rj^ TT- 

if[qi,t)...ip{qp+2,t) 

Then we obtain the following system of ODEs, 
at p 

(10) 

+ Hip,q,,q,)'f{q,t)v<P+'Hci,t), 

where q = (gi, . . . ,(?p+2)- 

Now consider the set flz defined by, 

Qz = \ w eC^ : w{0)> Cp,Amax{n2 \Re {w {n))\ |/m(u;(n))|} 

' n£Z 

To prove the lemma we must show that if i/jz = [i'z (n) ] belongs to ftz, so 

\ / nez 

does 

v{-,t) = {ij (ri,t))„g2, 

the solution to the ODE system pUj) with initial condition ipz, as long as it is 
defined {v{-,t), defined by the ODE system ((TU| . is a trajectory in C^, and what 
we want to show is that once a trajectory enters flz, it never leaves). In order to 
do so, we must show that whenever v belongs to flz up to time t = t, then 

dv , . f dv , . 



dt \dt J 

points towards the interior of ^z (in this case, r can very well be 0, and then we 
have V — ipz, which by hypothesis belongs to ilz)- 

Proving that ^ (r) points towards the interior of D,z whenever v (r) belongs 
to its boundary is a consequence of the fact that for Cp^x conveniently chosen the 
following inequalities hold (see Section 2 in [5]) 

(11) > 

J2qeAonZp+^ H {P, 91, 92) * (<?!, • • • , qp+2, t) ((ji, . . . , qp+2,r) 
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(12) > 

Eqe>t„n2p+2 {p, 92) ^ • • • , '?p+2, T") (qi, . . . , qp+2, t) 

So, in what foUows we will show that for a good choice of Cp^\, inequality (|12p 
holds whenever v belongs to the boundary of flz . The same reasoning can then be 
applied to prove inequality (jlip under the same circumstances, which would then 
prove the lemma. 

Now, if V belongs to the boundary of flz, then it holds that 



Cp,A max \\nf \Re{v {n,T))\ ,\nf \Im,{v {n,T))\\ < [ v{9,T)dd, 



^\nr \Ke{v{n,T))\Anr \lm(v{n,T))\)- < ^ 

and that there is an n e Z such that 

Cp,A ||w|l2 = Cp^A l^e {v (n, r))| = u (0, r) 

or the same, but for the imaginary part. Under these assumptions, if we write 
M = II I'll 2, the righthand side of inequality (IT^ is bounded above by 

/(cp,a)A2mp+2, 

where / is a polynomial of degree at most p, whereas the lefthand side is bounded 
from below (in absolute value) by 

P+l ( ;^2 _ P+2 7 \ ^p+2 

[ pn^ k{0,tr+'\n\ ' 



and hence, as 



pn^ k{0,t)P+^\n 



> 



(which can be achieved as long as A > y^^: choosing 7 > small enough), by 
taking Cp,\ > large enough, the lemma follows. The reader should have noticed 
also, that a judicious choice of Cp,\ > implies that k (0, t) is increasing (which is 
a consequence of inequality (fTTjl ). so the choice 7 > only depends on ip (0). 

□ 



Remark 2.1. The method of proof of Lemma \2.1\ gives a way to estimate Cp^x. 
Indeed, if p — 1 (the case of the curve shortening flow), we can take 

64A2 



Cl, A 



A2 



Although the Trapping Lemma is proven for finite dimensional approximations 
of the estimate given is strong enough so it "passes to the limit", i.e., it holds 
for solutions to ([5]), provided that the initial condition satisfies the hypothesis of 
Theorem 1 2 . II - the details are left to the reader. Hence, from now on our estimates 
are given for solutions to ([2]), and in consequence we drop the dependence on Z. 

Notice that the Trapping Lemma implies that solutions to ([2]) for initial condi- 
tions that satisfy the hypothesis of Theorem 12.11 are analytic in space. Also, from 
the Trapping Lemma, we can conclude the following useful estimate: for a solution 
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k to ^ with initial condition ip, which satisfies the hypothesis of Theorem I2.H we 
can find constants C, /i > such that the estimate 

T 



k (n, t) 



< Ce^^l"!, for t> ^, 



holds. Obviously C and fi may depend on the initial condition. The interested 
reader can compare this result with the work of Ferrari and Titi in 6] , where they 
prove analiticity results for certain semilinear parabolic equations in the d-torus. 

Also, from the Trapping Lemma and the ODE satisfied by k (0, t) we obtain the 
following result on the blow-up behavior of k (0, t). 

Lemma 2.2 (Blow-up Lemma). Let T > be the blow-up time of a solution to 
(0). Under the hypothesis of Theorem \ 2.1[ for every 77 > there exists a to > 
such that 

1 1 
fc(0,t)>' ^ ^ ^ " 

for all te {to,T). 



P+IJ {T-t)~^ 



Proof. Using the equation satisfied by k{0,t), by the Trapping Lemma it can be 
shown that k{0,t) satisfies the differential inequality 

4fc (0, t) < -k (0, tf+^ + Ak (0, tf , 
at p 

where A is a constant. Notice also that the Trapping Lemma implies that fc(0,t) 
blows up: since k blows up, and k {n, t) is conveniently bounded for n ^ 0, fc (0, t) 
must blow up. Hence, for every 77 > 0, there is a to > such that 

j-fc(0,0 < , ^ , -fc(0,t)P+^ forah to>0. 
dt {^-V)P 

The Lemma follows from integrating this differential inequality. 

□ 

Now we can give a first estimate on the rate of decay of the Fourier wave numbers 
of solutions to ([2]). 

Lemma 2.3. There exists cq > which depends on A such that if t > ^ > then 
there is a constant constant b > such that for any < e < eo, for n 7^ 0, the 
following estimate holds, 

T 

k{n,t) < 6e"^l"l (T-t)' whenever t>-. 

Proof. Notice that a solution to the infinite dimensional ODE system in Fourier 
space corresponding to equation ^ (whose finite dimensional approximations are 
described by (|8])) can be written as, 

(13) Hn,t) = fc(n,r)e-(^^"'-^)/*^("^^)''^^^^ 



+ 



J2 Hip,qi,q,)k<P+^Hq^,...,qp+2,t) ds. 
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From the Trapping Lemma, we can estimate the nonUnear term in the previous 
expression as 



J2 Hip,q,,q2)k*^P+'Hqi,---,qp+2,s) 
q.eA„ 



By the Blow-up Lemma, given 77 > there is a (5 > such that 'd t > T — d > 



then 



k {n, t) 



< 



k{n,T -5) 
+ {T-t) 



+ 



1 



T-5 {T-S 



C?S, 



where 



a (A, n,p) = A n 



2 2 P+2\ P 



P J P+^ 



With p fixed, by taking 77 = i, since a (A, n,p) > 0, we obtain a bound, 

<c(^\k{n,T~6)\ (T - t)' + e^^l"! [T - t)^^ , 
for any < e < min {^ct {X,n,p) ,5}, which proves the lemma. 



□ 



The previous lemma already implies our stability results. However, in order to 
obtain sharp estimates on the rates of uniformization of solutions to ([2]), and to 
finish the proof of Theorem l2.H Lemma [2.2l does not suffice. In fact, we can improve 
a bit on Lemma [221 So we have: 



Lemma 2.4. There exists a to > such that for all t G {to,T) the following 
estimate holds 



k{0,t) > 



P 



p+1 



1 



(T-t) + {T-t)^+^^'^' 



Proof. As in the proof of Lemma 12.21 the following differential inequality holds 

4fc(0,i) < -k{0,tf+^ + Ak{0,tf , 
at p 



for a constant A > independent of t. So using Lemma [2.21 we obtain the differ- 
ential inequality 

^—^^k{0,t) < - + C{T-t)^^ , 

kiO,tf+^dt ^ ' ' - p ^ ' 

which after integration gives the desired inequality. 

□ 



We are ready to improve the estimate on the decay of the Fourier coefficients of 
solutions to ([2]), i.e., the estimate provided by Lemma [^751 In order to proceed, we 
use the previous lemma to estimate the integral 



p+1 
P 



k{0,t) 



P+1 
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from below. Indeed, from Lemma 12.41 since t and t are close to T, using Taylor's 
Theorem, a simple computation shows that 



/ > 



■ ds — — los 



T-t 

T-T 



T - s + (T - s, 

and using this and (IT^ to estimate k {n, t) from above, yields 

'ji _ ^\ a(A,n,p) 



0(1), 



(14) 



k (n, t) 



< C 



k{n,T-5) 



t 

T-S 



T- 



a{X,n.p) 



X! \h {p,qi,q2)k*'^P+^^ (gi,...,qp+2,s) 
qe-4„ 



ds. 



where 



a(A, n,p) —[An 

V P 

Let us now improve on the estimate given by Lemma 
from Lemma 12.31 into (IT4|) , we get 

a( A,n,p) 



P 



P+1 

If we introduce the bound 



k (n, t) 



< 



k{n,T-S) 



T-t 



(T — s f'^ ds, 



Its [T - s)"(^'"'f ) 
with < /i' < /i, and from which we obtain the estimate 

k {n, t) < C'e-^'l"l (T - i)-i"{"(^."^P).i+2e} . 

Using this new bound and plugging it into (|14p. we improve again our estimate on 
k {n, t): 

k (n, t) < C"e-^"l"l (T - i)™{"(^.».p).3+4.} (0 < ^» < 

Finally, it should be clear that if we repeat this procedure a finite number of times, 
we arrive at 



(15) 



kin,t) < i?e-''l"l (7^ - "7^0, 



which in turn implies estimate (|4]), i.e., the conclusion of Theorem 12. II 

3. A FEW COMMENTS ON THE STABILIZATION RATE TOWARDS THE CONSTANT 
STEADY STATE SOLUTION OF THE NORMALIZED p-CURVE SHORTENING FLOW 

In this section we study the exponential stability of the steady solution u = 1 
of the normalized p-curve shortening flow. Recall that the normalized p-curve 
shortening flow is obtained from the unnormalized p-curve shortening flow by the 
process described in section [2l To be more precise, we have that the normalized 
p-curve shortening flow is equivalent to the Boundary Value Problem 



du 



(16) 



^ ^ pu^^p ^ + pu'^+p - pu in [0, ^] X (0, cx)) 
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with periodic boundary conditions. Recall that this normalized equation is obtained 
from ([2|) by the rescaling and change of time parameter given by 



u{9,t) 



P 



p+1 



{T -t)^ k{0,t) , 



p+1 



log 1 



T 



The reader should have noticed that Corollary 12.11 does not give a rate of conver- 
gence of the solution to the normalized flow towards the steady solution u = I. In 
order to provide rates of convergence towards the steady solution, we will analyze 
the behavior of k (0, t) in the case of the unnormalized equation. First observe that 
the following lemma holds, and that its proof is an obvious modification of the 
proof given in Lemma 12.41 Again, we are under the hypothesis of Theorem 12.11 

Lemma 3.1. There is a to > such that iftCz {to,T), then we have the estimate 



kiO,t) < 



p+1 



1 



[T -t)- (T-<)^+5^ 



If we let 



u{Q,t) 



P 



and then compute 

u(0,i) - 1 < 



p+1 



(T-i)— fc(0,i), 



(T-t)- (T-t)^^^ 



[T -t)- [T -tf^^ 



factoring out from the numerator (T — and applying Taylor's Theorem we 
obtain 



[T-ty^ 



— {T~t)~ <C{T ~ t)~ 



M(0,i) - 1 < 

{T-t) - [T -tf^^ 
In the same way, by using Lemma 12.41 we obtain an estimate 

l-u{Q,t)< C{T-t)^ . 
This implies that for the normalized flow ([T6l) holds that 
A 



2^ Jo 



u [9, t) dO-l 



< Cexp(-2r). 



Using this estimate and Corollarv l2.1[ we obtain the estimate 
\\u{0,t) - l||^oo([o^2,]) < Cexp(-WT), 

u! = min |A^p — p — 1,2} . 
A^p - p - 1 < 2 



(17) 
where 

Notice that 

whenever 



p + 2 



< A < 
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SO it is in this case that we obtain a rate of decay towards the steady state cor- 
responding to the first eigenvalue of the elliptic part of the normalized p-curve 
shortening flow. However, if we use estimate (jf 5p in the proof of Lemma 12.41 we 
can bound 



with 



<C{T-t) 



k{0,tf 



P 



P 



1' 



then from the ODE satisfied by k{0,t), we can obtain the improved bound from 
below 

p+i I 



k{0,t) > 



P 



p+1 



1 7 
P+1 



and also the corresponding bound from above. Hence, proceeding as we just did, 
we arrive at an estimate 



— / uie,T)de-i 



< Cexp(- {2X^p~p) t) 



from which follows that (ITTI) holds now for all A > with 

cj = min {X^p — p — 1, 2X?p — p} = X^p — p — 1. 
So we finish with the following result. 

Proposition 3.1. Let X > ^ let iJj > and assume that it satisfies the 

hypothesis of Theorem \2.1[ Then the solution u{9,t) to I116\) corresponding to the 
solution to (0i with tjj as initial data, satisfies an estimate 

\\u{e,T) - l||i=„[o,^] < C'exp(- (A^p-p- 1)t) , 

where C > is a constant that depends on tjj, p and X. 
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